Dynamical domain walls and spin-Peierls order in doped antiferromagnets: 
evidence from exact diagonalization of small clusters 
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The hole-doped 2D t-J model is studied by exact diagonalization on a 5 x 4 cluster which, unlike 
the standard tilted square clusters, can in principle accomodate an antiphase domain wall. For hole 
concentration 10% and J/t > 0.5 the ground state energy/site is lower than the conventional tilted 
square V20 x v^-cluster. In the ground state two holes form a loosely bound pair pinned to an 
antiphase domain wall. The dynamical density correlation function shows sharp quasiparticle-like 
peaks, reminiscent of the 'holons' in ID chains, which suggest the existence of soliton-like propagat- 
ing domain walls. The dynamical correlation function of the bond-singlet operator has a low-energy 
peak structure characteristic of columnar spin-Peierls order, the dynamical spin correlation function 
shows an intense and isolated 'resonance peak' near (7r,7r). 
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presence of stripe-like structures [3| which were 
pf^icted early on by mean-field studies of Hubbard-like 
rd^elsi^, is the basic idea underlying much of the re- 
C(MiJ theoretical work|l H S E 

on cuprate supercon- 
drartors. An interesting question is whether the 'pure' 
2^t-J model, which still represents the simplest the- 
oijalical description of the Cu02 planes, develops any 
stQje-like structures, either in its ground state or in low- 
qiierg y excited states. The most conclusive evidence for 
the spontaneous formation of such stripes in the 2D t- 
JrTTiodel is probably the density-matrix renormalization 
geetip (DMRG) study of White^nd Scalapino 01 and the 
r^lkted work by Martins et al. 9] . Exact diagonalization 
) of small clusters on the other hand, so far has 
produced really compelling evidence for hole stripes, 
ast in the physical parameter regime. Evidence for 
S(titpe-like hole-density correlations from ED calculations 
lS&>been reported by Prelovshek and Xotos[Tlj|. but only 
fq^elatively large values of J/t > 1.5. This is astonish- 
ing because all effective interactions in the t-J model are 
icted to be short ranged, whence any strong tendency 
rm stripe-like structures should make itself felt even 
lall lattices - unless it is suppressed by finite-size ef- 

A^lready stressed by White and Scalapinoj^ a potential 
s(i^ce of such finite-size effects are the boundary con- 
ditions imposed by the cluster geometry. Most studies 
t^^ate have been performed on tilted-square clusters 
wWi periodic boundary conditions (PBC) , which are cho- 
seir such as to accomodate the Neel-type order parame- 
ter in the undoped system. On the other hand, such a 
geometry is precisely the wrong one when an antiphase 
domain wall is present in the ground state. The mere 
geometry of the cluster then can frustrate the domain 
wall, thus artificially enforcing a homogeneous ground 
state. An ideal system to check this would be the 5x4- 
cluster with PBC, which on one hand is appropriate to 
accomodate an antiphase domain wall parallel to the y- 
direction and on the other hand can be compared di- 




FIG. 1: Difference AEo between the GS energy of the 5x4 
and the \/20 x \/20 cluster, as function of J/t. 



rectly to the tilted square 20-site cluster. As will be 
demonstrated, this cluster shows very clear and unam- 
bigous evidence for a stripe-like domain wall not only in 
the ground state, but also in the form of excited states 
corresponding to a soliton-likc propagating domain wall. 
The 2D t-J model thus has a very strong intrinsic ten- 
dency to form stripes - if it is allowed to do so. Moreover, 
the formation of stripes seems intimately related to spin- 
Peierls-like columnar singlet order P. IT2|. 
The Hamiltonian of the t-J model reads 



H = -t 



H.C 



niUj 



Here denotes summation over nearest neighbor 

pairs, Ci^a- — Ci^cr{l — ni,g) and Si and rii denotes the 
operators of electron spin and density at site i, respec- 
tively. We begin with a comparison of the ground state 
(GS) energy for the 5x4 cluster and the tilted square 
20-site cluster. One might expect that this can give an 
indication, as to whether a given cluster has the proper 
geometry to describe the system. Inappropriate bound- 
ary conditions introduce 'frustration' thereby raising the 
energy. For example, at half filling the GS energy in 5 x 4 
is —1.165 J/site, whereas in the square-shaped cluster it 
is — 1.191J/site, the higher energy for 5 x 4 obviously 



2 



being due to the frustration of the (quasi-) Neel order 
along the odd-numbered side of this rectangular cluster. 
Interestingly enough, this changes when holes are added. 
Figure n compares the GS energy of the two clusters with 
2 holes as function of J/t. For J/t > 0.5 the rectangular 
cluster indeed has the lower energy, indicating that for 
the doped case the different boundary conditions are at 
least as reasonable to represent the bulk system as is the 
square cluster. All data to be presented below have been 
obtained for the value J/t — 0.5 and 2 holes. Table U 
shows the static density and spin correlation functions 
for the 5x4 cluster with 2 holes. These are defined 
as the GS expectation values f/£)(R) = X^^i'^i^j+R.) 
and 55 (R) = 'h'^j^^i ' ^j+R-)- There is pronounced 
anisotropy in (?d(R) , particularly so at short distances, 
which immediately suggests a 'hole stripe' in y-direction. 
The fact that the hole density correlation is stronly direc- 
tional might seem to suggest that this state is related to 
the p-like pairing states, which usually form the first ex- 
ited state above the dx^-y2-\\ke ground state (and some- 
times in fact the ground state itself) ^^of two holes in 
square clusters. These p-like pairs, however, are spin- 
triplets whereas the present ground state is a spin singlet 
which moreover is even under reflection by both the x- 
and y-axis. This state therefore is qualitatively very dif- 
ferent from those seen so far in square-shaped clusters. 
Pronounced anisotropy is also seen in the spin-correlation 
function, where in particular the nearest neighbor spin 
correlation in y-direction (i.e. parallel to the stripe) 
exceeds the one in x-direction by 50%. While some 
anisotropy in the correlation functions is to be expected 
solely due to the rectangular shape of the 5x4 cluster, 
the edges differ in length by only 20% which seems to be 
rather small to explain the strong anisotropy. 
To further clarify the nature of the ground state, Figure 
121 shows the 'representatives' of those basis states, which 
have the largest weight in the GS wave function, labelled 
by their coefficients in the GS. Taken together, these 9 
states exhaust 14.1% of the GS wave function - which is 
a lot, given that there are all in all 59000 basis states for 
this system. For clarity we mention that a basis state 
for the ED scheme is created from the 'representative' by 
a) translating it by all N different lattice vectors of the 
cluster, b) applying all iV' point group operations of the 
little group of the total momentum k c) applying the time 
reversal operator and d) adding up the resulting 2NN' 
states with the proper phase factors so as to create a state 
with prescribed momentum, point group symmetry and 
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FIG. 2: Hole-spin configurations of the representatives with 
the largest weight in the GS wave function. Dashed and solid 
boxes indicate the hops and spin flips to created the state 
from the 'seed state' in the top left panel. 



parity under time reversal. All states in Figure [3 show 
a domain wall separating two perfect Neel states. In the 
configuration with the largest weight (top left) the two 
holes form a pair parallel to the domain wall. All subse- 
quent states can be generated from this 'seed state' either 
by hole hopping, by a quantum spin flip along the domain 
wall, or by a combination of both. This is very much rem- 
iniscent of 'string picture' theories for ho les p opulating a 
domain wall between two Neel states [T^ hj and m par- 
ticular highlights the importance of charge fluctuations 
transverse to the stripe 1^, ■ Taken together the data 
presented so far show that the 5x4 cluster has a quali- 
tatively new (as compared to the tilted square clusters) 
ground state: a loosely bound hole pair pinned to an an- 
tiphase domain wall of the staggered magnetization, with 
strong quantum spin fluctuations along the fault line. 
Having established the nature of the ground state, we 
consider the excitation spectra of the system. As a first 
step, we address the existence of dynamical domain walls. 
Here the appropriate quantity to look at is the dynamical 
density correlation function (DCF), defined as 



D{q,u;) 



1 



(1) 



where «.„ 



'7m nje^"^'^^ is the Fourier transform 



TABLE L Static correlation functions for the 4x5 cluster. 



"1 ~ Vn 

of the electron-density operator rij . Figure |3| compares 
D{q,oj) for the 5x4 and the tilted square 20-site clus- 
ter. Whereas the DCF for the tilted square cluster is 
essentially incoherent . the DCF for the rectangular 
cluster shows intense and well defined peaks at the lower 
edge of the spectra, particularly so at q = (27r/5, 0) and 
q — (47r/5,0). In fact, it is tempting to compare the 
DCF for the 4x5 cluster to the one for ID chains. Fig- 
ureOalso shows the DCF for an 11-site t-J chain with one 
hole. The DCF here really is an almost ideal single-peak 
spectrum, with a dispersion that can be fitted very accu- 
rately by the expression e{q) = 2t{l + cos{q)). Obviously 
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FIG. 3: Left: DCF for the 4 x 5 and V20 x V20 cluster. The 
upper/lower momentum label refer to 4 x 5/\/20 x \/20. 
Right: DCF for a single hole in a 1-D 11-site ring at J/t = 0.5. 



FIG. 4: Left: Ba{q_,uj) and 5'(q, cj) in the 4x5 cluster 
with two holes. Right: The same correlation functions for the 
\/20 X \/20 cluster, ^(q, lj) is multiplied by 0.5. 



these peaks originate from the particle-hole excitations 
of the single 'holon' in the system. Since the holon is 
nothing but a domain wall, this suggests that the sharp 
peaks in the 5x4 DCF for momenta (g^, 0) also originate 
from the propagation of the domain wall as a whole. The 
stronger q-dependence of the peak weight in the 2D clus- 
ter is probably due to the finite width of the domain wall: 
whereas the holon really is a point-like object, making its 
density structure factor gjj (q) flat in g-space, the domain 
wall in 2D has some extension in x-direction, whence its 
structure factor must have a g^j-dependence. 
Next, we focus on the anisotropic spin correlations. 
Read and SachdevjlTj have proposed a spontaneous spin- 
Peierls-like dimerization to occur as a general feature of 
S — ^ antiferromagnets in 2D - which would explain the 
anisotropy of 35 (R). The appropriate quantity to check 
this hypothesis is Ba (q, , the dynamical correlation 
function (defined as in (Q) of the bond-singlet operator 
Ba,^ = (S, • S,+e„ - ^^^) e^"^-^^. Thereby 

a G [x, y] denotes the direction of the singlet-bond. Also 
of interest is S'(q, cj), the dynamical correlation function 
of the spin operator 5q = J2j Sje^'^'^' . These spec- 
tra are shown in Figure 21 To begin with, Ba{q^,uj) in 
5x4 has a prominent low energy peak for both sin- 
glet directions, a — x,y, at q = (27r/5,0). This peak 
(as well as the similar peaks at (0, ^) and (47r/5,0)) is 
simply a replica of the intense low energy peak seen at 
the same q in the DCF. Just as the density operator 
itself, Ba,q is a number-conserving spin singlet, whence 
Ba (q, uj) probes exactly the same final state manifold 
as the DCF. Since the bond singlet operator (partially) 
couples to fluctuations of the electron density, Ba{q,io) 
must pick up the signal from the propagating domain 
wall. This interpretation is confirmed by the fact, that 
the excitation energies E^, — Eq of the respective peaks 



in the DCF and Ba (q, oj) agree to computer accuracy, 
i.e. these peaks originate from the same final state 
More interesting therefore is the second high-intensity 
and low energy peak, at (0,7r). This peak appears only 
in -Bj,(q, to), that means for bonds parallel to the domain 
wall. The wave vector (0, tt) for bonds in y-direction cor- 
responds exactly to the Spin-Peierls order found by Read 
and Sachdevp^Tj and proposed to be a general feature of 
doped antiferromagnets by various workers!^ [ij. This 
peak has a lower excitation energy than the dominant 
peak in the spin correlation function S{q,u;) at (^,7r), 
indicating that spin-Peierls order is the most likely in- 
stability of the system. S{q, uj) itself is quite remarkable, 
especially when we compare it to the tilted square clus- 
ter. There S'(q, w) shows a series of low energy peaks 
near q = (tt, tt), which suggest a relatively smooth dis- 
persion with a shallow minimum at the incommensurate 
wave vector (^, ^). S{q,Lu) in 5 x 4, on the other hand, 
consists of a rather diffuse high-energy continuum and an 
isolated low-energy peak at (^, tt). The dispersion of the 
spin excitations (if it exists) must be rather discontinuous 
near this momentum. Such a spin-excitation spectrum is 
very much reminiscent of the 'resonance peaks' observed 
in neutron scattering experiments on the superconduct- 
ing state 10. In the framework of the spin-Peierls sce- 
nario, the peak should be interpreted as follows: start- 
ing from a columnar spin-Peierls state with singlets in 
y-direction, a Neel-ordered state can be generated by con- 
densation of (Bosonic) bond-triplet excitations [13 with 
momentum Q — (7r,0). In the present case of an an- 
tiphase domain wall, this should be replaced by Q = 
(^, 0). The spin operator ^-j then creates and anni- 
hilates the condensed triplet-excitationsjl^, so that the 
intense low-enery peak at (^,7r) is simply the finger- 
print of the (quasi)-condensate of triplets, which creates 
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FIG. 5: Mechanism for allignmcnt of hole pairs. 

antiphase Neel-like spin correlations out of the columnar 
spin-Peierls state. To conclude our discussion of the spin 
dynamics, we note that By{q, uj) for the \/20 x \/20 clus- 
ter also does show a low-energy peak at (0, tt) - the ten- 
dency towards spin-Peierls instability thus really seems 
to be quite generalj^ 0, 0|. Since the geometry of 
square cluster is unfavourable for an antiphase domain 
wall, however, this does not make itself felt as clearly as 
in the 5x4 cluster. 

In summary, it has been shown that a rectangular clus- 
ter of the t-J model whose geometry does not explicitely 
frustrate an antiphase domain wall, shows rather differ- 
ent behavior than the tilted square clusters used conven- 
tionally for exact diagonalization. The GS in this cluster 
shows the clear and unambiguous signatures of a stripe- 
like domain wall, a comparison of the GS energies shows, 
that such rectangular clusters are at least as well suited 
to describe the bulk system, as are the square clusters. 
Clearly, this does not prove the existence of stripes in 
the thermodynamic limit - just as a square-shaped clus- 
ter will tend to suppress a domain wall, the rectangular 
one will favour it - but given the previous evidence for 
stripes found by White and Scalapinoj^ on much larger 
systems, it is quite plausible that the stripes in the rect- 



angular cluster have a similar structure. 
Stripe-like structures in the rectangular cluster then ap- 
pear not only in the ground state - rather, there is evi- 
dence for excited states where apparently the finite total 
momentum of the state is carried by a soliton-like prop- 
agating domain wall - which is very much reminiscent 
of the 'meandering domain wall' scenario put forward 
by Zaanen et aZ.,^]. The data also show the presence 
of quasistatic columnar singlet order |^ ^| in the sys- 
tem, which seems to be either a prerequisite for or an 
immediate consequence of the stripe formation. Taken 
together, the above results suggest the following scenario 
for doped antiferromagnets: the spontaneous breaking of 
the point group symmetry of the lattice^3| by formation 
of bond-singlet order produces ladder-like patterns in the 
spin background of the system, which then serve as 'race- 
tracks' for hole-pair-like solitons. Condensation of bond- 
triplets with a condensation amplitude that changes sign 
across a soliton introduces strong antiferromagnetic cor- 
relations with opposite staggered magnetization on the 
two sides of the soliton. Under these circumstances, it be- 
comes energetically favorable for the solitons on different 
ladders to form a line, because then the staggered mag- 
netizations on the different 'ladders' can allign. Solitons 
which propagate away from the stripe create a track of 
magnetic frustration (see Figure|SJ), much as a single hole 
in an antiferromagnet, which creates an effective poten- 
tial that alligns the pairs. In this way, one arrives at the 
picture of a fluctuating domain wall^, made of loosely 
bound hole pairs in a background provided by the colum- 
nar singlets. An alternative point of view, which also fits 
the data very well (see especially Figure would be the 
accumulation of hole pairs along the fault line separating 
two Neel states with oppposite staggered magnetization, 
as proposed early on by White and Scalapino[2H. 
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